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^  This  paper  is  the  author's  contribution  to  the  volume  *Leonardt  Euler- 
Gedenkband  1983*  to  be  published  in  Basel,  Switzerland,  in  1983  in  honor  of 
Euler's  bicentenial.  It  is  mainly  devoted  to  the  exponential  Euler  spline 
SJ^(x/t)  of  degree  n  to  the  base  t,  and  also  sketches  in  §7  their  role  in 
cardinal  spline  interpolation.  It  also  presents  two  new  items:  1.  In  §3  a 
simplified  derivation  of  the  recursive  relation 

Sn(x,-t)  =  I**1  S^U/Udu  //’  Sr_1  (u  ,”t)du  (n  >  2) 
is  given,  a  relation  already  discussed  in  Reference  [10).  2.  In  §6  the 

approximation  of  the  exponential  function  2X  by  Sn(x/t)  is  made  more 
effective  by  a  preliminary  subdivision  of  the  interval  [0,1]  into  2r 
parts.  By  this  device  our  approximation  becomes  competitive  with  the  modern 
approximations  of  2X  in  [0,1]  by  rational  functions.  The  paper  has  two 
aims:  1 .  As  a  tribute  to  Euler,  2.  To  make  the  exponential  Euler  splines 

S£(X;t)  better  known. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  is  the  author's  contribution  to  the  volume  "Leonardt  Euler  - 
Gedenkband  1983"  to  be  published  in  Basel,  Switzerland,  in  1983  in  honor  of 
Euler's  bicentenial.  It  has  two  aims:  1.  As  a  tribute  to  Euler,  2.  To  point 
out  Euler's  contributions  to  the  subjects  of  the  title. 


i. 


,  *  ;  on_. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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EULER'S  CONTRIBUTION  TO  CARDINAL  SPLINE  INTERPOLATION: 

THE  EXPONENTIAL  EULER  SPLINES 

X.  J.  Schoenberg 

Introduction*  In  my  monograph  [8]  of  1973,  dedicated  to  Euler,  I  already  discussed 
the  subjects  of  the  title.  On  the  occasion  of  the  bicentenial  of  Leonardt  Euler  we  present 
here  an  outline  of  these  results,  which  seem  to  fit  well  in  what  we  think  of  as  Eulerian 
Mathematics. 

Our  main  subject  are  the  exponential  Euler  splines.  In  §1  we  define  them,  and  §2 
shows  their  close  connection  with  the  Eulerian  polynomials.  In  §3  we  derive  in  a  simpler 
way  a  recursive  construction  already  described  in  [10].  §§4  and  5  show  that  the 
exponential  Euler  splines  of  base  t  converge  to  the  exponential  function  t*  as  their 
degree  tends  to  infinity.  §6  presents  an  application  to  the  computation  of  f(x)  *  2X. 
Finally,  in  §7  we  sketch  the  role  of  the  exponential  splines  in  the  problem  of  cardinal 
spline  interpolation. 

1 .  The=exgonential=t|uler_ splines.  We  need  a  few  definitions.  Let  *  (s^ ( x )  } 

denote  the  class  of  cardinal  splines  S^t*)  °f  degree  n  (>  1).  This  means  that  Sn(x) 
reduces  to  a  polynomial  of  degree  <  n  in  each  unit  interval  (v,  v+1)  (V  €  Z),  with  the 
strong  restriction  that 

(1.1)  S  (x)  e  Cn_1(R)  . 

n 

In  particular  S^x)  e  Sj  means  that  S1(x)  is  a  continuous  piecewise  linear  function 
with  possible  vertices  (or  "knots")  at  the  integers.  Early  in  this  century  it  was  found 
convenient  to  represent  S^x)  as  a  linear  combination  of  shifted  versions  of  the  "roof- 
function" 

Q2(x)  =  x  in  (0,11,  *  2  -  x  in  (1,2),  =  0  elsewhere  , 

so  that 

00 

s t(x)  *  I  Cjc2(*  “ 

represents  uniquely  every  element  of  Sj 
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(1.2) 


This  extends  to  the  class  $  in  terms  of  the  forward  9-spline 


n¥  1 


*n+r 


1*1  -t  i  . 


V-0 


where  u+  -  max(0,u).  Like  Q2<x)#  ***•  B-spline  Qn+1(x)  hee  remarkable  properties: 

(1.3)  Qn+1(x)  >  0  in  (0#i*4l)#  -  0  outside  • 

Moreover,  it  is  bell-shaped  in  (0,n+1)  and  symmetric  in  its  midpoint#  i.e* 

(1.4)  Qn+t<n  +  1  *  x>  *  • 

clearly  ©^(x  “  j)  e  Sn  for  all  integers  j#  and  these  are  the  elements  of  SR  of 
least  support.  Again#  every  element  S(x)  6  $n  a<taits  a  unique  "standard"  representation 
of  the  form 

<1.5)  S<x)  -  l  e  O  («  -  j)  . 

_•  J  n+1 


Definition  1 ;  The  exponential .splines  og,base  t.  Let  t  be  a  real  or  complex 
number  +  0,  and  let 

m 

(1.6)  *n<x,t>  "  I  ^  • 


We  call  this  function  the  exponential  spline  of  degree  n  and  base  t. 
Clearly 

♦  ( x*  1 1 1 )  ■  £  t^Q  - (x+i- j )  -  Z  t^Q  .t(x-j)  -  t  •  (xit)  • 

n  n-M  nti  n 

Using  the  representation  (1.5)  and  its  unicity#  it  is  essily  shown  ( (7#  Lemma 

most  general  solution  of  the  functional  equation 

(1.7)  S(x  +  1)  -  t  8(x) ,  wh.r.  S(x)  «  Sn 


is  given  by 


(1.8) 


S(x)  -  c*d  (xit),  (C  is  a  constant)  . 
n 


If 


(1.9) 

let  us  try  to  interpolate 

(1.10) 


|t|.ia,  -w  <  a  <  «,  t  t  0 


the  exponential  function 
^x 


tix  .ia  x 


at  the  integers  by  the  function  (1.8)  so  that 

S(v)  -  tV  for  all  Integers  v  • 


2) )  that  the 


(1.11) 


-.J 


Because  of  (1.7)  it  suffices  to  determine  the  constant  C  in  (1.8)  so  that  S(0)  «  1.  The 
answer  is  clearly 


(1.12) 


S(x) 


♦  (x; t ) 
n 

♦  (0»t)  ' 

n 


but  this  is  possible  if  and  only  if  ♦  (Ojt)  J*  0.  When  this  holds  is  easily  decided,  for 

n 

by  (1.6)  and  (1.4)  we  have 

*n(0jt)  -  E  t3Cn+1(-j)  =■  E  t3fin+1(n+1+j)  , 

and  setting  n+j  -  v  we  find  that 

(1.13)  *n(0,t)  =  l  tV'nQn+1(v+1)  -  t'n  l  Qn+1 (v+1 )tV  . 

The  result:  The  interpolation  (1.11)  with  S(x)  of  the  form  (1.8)  is  possible  if  and  only 
if 

n-1  . 

(1.14)  n  (t)  =  nl  l  Q  ( j+1 )tJ  t  0  . 

n  „  n+1 


The  polynomial  n^(t)  defined  by  (1.14)  is  called  the  Euler-Frobenius  polynomial.  It  is  a 
reciprocal  monic  polynomial  having  integer  coefficients  and  having  only  negative  and  simple 
zeros 

,,.15)  An_,  <  An2  <..•<  X2  <  A,  «  0,  . 

Definition  2:  The_exponential  guIgr_ splines  SR(  x; t ) .  Assuming  that  II^( t )  ^  0, 
hence  that 


(1.16) 
we  define 
(1.17) 


♦  (Oft)  t  0  , 


♦  (Xft) 

sn<x,t)  “  TTo 7T7 


To  summarize  Sn(x;t)  is  the  unique  cardinal  spline  interpolant  of  the  exponential 
tx  satisfying  the  functional  equation 


(1.18) 


Sn<x+1*t)  -  t  Sn(x;t>,  (x  e  R) 


-3- 


2.  The  construction  of  S  (xjt)  in  terms  of  Eulerian  polynomials.  How  do  we 
construct  Sn(xjt)?  Clearly,  its  expression  by  (1.17)  is  too  laborious.  This  is  where 
Euler  comes  in.  Following  Euler,  we  define  the  an(t)  by  the  expansion 

•  a  ( t ) 


(2.1) 


t-1 


t-e*  0  n 


c  n  n 
l  nl  z 


The  an(t)  are  rational  functions  of  the  form 


(2.2) 


a  ( t) 


n  (t) 


u-ir 


where  II  (t)  are  the  polynomials  (1.14).  For  a  proof  see  (7,  Lemma  7  on  page  391 J.  The 
n 

II  (t)  may  also  be  defined  by  Euler *8  expansions 
n 


(2.3) 

He  find  that 


n  (t) 

n 


(1-t)n+1  VO 


l  <v+1)V  . 


nQ(t)  -  n^t)  -  i,  n2<t)  -  t  ♦  i,  n3<t)  -  t  +  4t  +  i  , 

K.(t)  -  t3  ♦  lit2  ♦  lit  +  1,  n.(t)  -  t4  +  26t3  ♦  66t2  +  26t  ♦  i  . 

4  3 


On  multiplying  (2.1)  by  exz  we  obtain  Euler's  generating  function 
(2.4) 

of  the  exponential  Euler  polynomials 

n-1  .  rnv  ,^».n-2 


^  a  *  A  (Xjt) 

t-1  kz  v  n  n 


(2.5) 


A  (x#t)  -  xn  ♦  (”)a  (t)x0’1  +  (”)a  (t)x”"2  ♦•••♦«  (t)  , 

n  ii  z  2  n 


which  evidently  form  an  Appell  sequence  (see  12,  Chap.  VII,  178]).  L.  Carlitz  [1]  writes 
An(x;t)  -  Hn(x|t)  and  calls  them  Eulerian  polynomials.  See  also  [1)  for  extensive 
references. 

The  coefficients  an(t)  a*»it  a  recursive  computations  Multiplying  (2.1)  by  t-e* 
we  obtain 


t 


•  ta  (t)  •  1  ♦  (™)a  ( t)  +  •••  +  a  (t) 

r  n  n  r  1  1  _  n  ^n 

a  nl  *  "  A 

0  0 


nt 


f 
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3rt*s 


and  by  identifying  coefficients  of  zn  we  obtain 

(2.6)  1  +  +  (”]a  (t)  +  •••  +  a  (t)  *  t  a  (t),  (n 

ii  z  a  n  n 


which  show  that 

(2.7)  a  (t) 

n 


t^T  {1  +  +  "•  +  On-i(t,}'  " 


1,2, 


)  , 


1,2,...)  . 


Let  us  remember  that  we  wish  to  construct  Sn(x;t),  and  that  we  may  exclude  the 
trivial  case  when  t  *  1,  because  evidently  Sn(xj1)  «  1  for  all  x.  We  ash:  What  can 
we  say  about  the  function  F(x)  defined  by 

(2.8)  F(x)  *  A  (x;t)  if  0  <  x  <  1  , 

n  —  * 

and  satisfying 

(2.9)  F ( x+ 1 )  *  t  F(x)  for  all  x  ? 

We  claim  that 

(2.10)  F{ x)  e  Cn” 1 ( R)  . 

indeed,  from  (2.10),  and  using  (2.9)  and  (2.8),  we  obtain  by  differentiation  of  (2.9) 
and  setting  x  «  0,  that  we  must  have  that 

(2.11)  A*V,(1it)  -  t  A*V)(0>t),  (w  -  1 , . . . , n- 1 )  for  n  >  1  . 

n  n  = 

However,  these  relations,  together  with  the  fact  that  An(xjt)  is  monic,  are  known 
characteristic  properties  of  An(xjt),  which  are  derived  from  (2.4)  by  v 
differentiations  with  ' espect  to  x,  subsequently  setting  x  *  0  and  x  *  1  in  the 
result.  By  our  result  (1.8)  concerning  (1.7),  it  follows  that 

(2. 12)  F( x)  «  C  ♦  (x/t) 

n 

Assuming  (1.14),  this  proves  that  Sn(xjt)  -  An( x; t)/An( 0; t)  in  [0,1],  hence  that 

(2.13)  S(x,t)  -  (xn  +  (”)a  (t)xn  1  ♦  •••  +  a  (t)}/a  (t),  if  0  <  x  <  1  . 

n  ii  n  n  *  * 


Remarks.  1.  In  (1,  page  256,  (4.5)]  Carlitz  already  defined  the  cardinal  spline 
F(x)  satisfying  (2.8)  and  (2.9). 

2.  The  continuity  requirement  (2.11)  has  been  recently  stated  as  a  general  principle 
concerning  the  solutions  of  certain  functional  equations  in  the  paper  [3]. 

3*  Sn(x?t).  This  is  the 

subject  of  my  recent  paper  [10],  with  the  modification  that  there  the  sequence 


(3.1) 


8^ (xjt) 


8j(*»t) 


s2(*+  y  t) 


S4(x*  yt) 

84(ft) 


is  recursively  constructed.  We  must  t  to  be  noa-naoativa.  t  +  0.  For  the  B-spline 

(3.2)  Qn(x)  -  T^T^TY  Z  <-DW(“)  (X  -  V,-' 

we  easily  verify  by  integration  and  su— itioci  by  parts  the  relation 

(3.3)  /f*  ^<«>*>  -  * 

For  the  exponential  spline  (1.6)  this  implies  that 


r1 

9  X 


♦  (uit)du  -  l  t3  r*  Q  (u-j)du  -  It3  /^3  1  Q  1 

n  1  j  51  n  j  “  n 

-  Z  tlVi(st_J>1)  “  1 1 


(u)du 


iH-r 


and  therefore 

(3.4)  t^^uitjau  -  t  #o(x|t)  • 

For  the  exponential  Buler  spline  this  implies  our  first  proposition 
X*  It  t  is  not  negative#  then 

(3.5)  Sn<xit)  -  8^  (uiOdu/Zg  S^CuitJdu,  (n  -  2,3....)  • 

This  is  remarkable  recursive  constructions  Starting  from  the  linear  Buler  spline 
Sj(x/t),  (3.5)  recursively  furnishes  all  higher  degree  8n(xit).  Also  notice  that  (3.5) 
does  not  depend  on  t  explicitly. 

We  also  need  the  following  result  established  in  (•»  Lecutre  2,  §5] t 
IX.  If 

(3.6)  t  -  |t|.ia,  -«  <  a  <  «,  t  *  0,  t  t  1  , 

which  implies  that  negative  values  of  t  are  excluded,  then 

(3.7)  #  (xit)  +  0  for  all  real  x  • 

n 

The  reason  given  in  (0,  loc.cit.]  is  as  follows!  1.  If  t  >  0,  then  the  curve  of  the 


-6- 


—  -  •  - 


complex  plane 


(3.8)  r  S  z  -  ♦  (X|t)  (-•  <  x  <  +•) 

n 

ia  clearly  contained  within  the  positive  half  of  the  real  axis.  2.  If  in  (3.6)  we  have 
0  <  a  <  v,  say,  then  the  curve  (3.8)  spirals  convexly  about  the  origin,  never  assuming  the 
value  z  -  0.  This  is  shown  by  induction  in  n. 

4.  A^serleg^expanglon^of  Sn(x;t).  Again  we  assume  that  (3.6)  holds.  Already  in  (7, 
§7]  we  derived  by  means  of  residue  theory  the  following  proposition 
III.  Let 

(4.1)  Y  =  log  t  *  log  |t|  ♦  ia,  hence  t  *  e^  , 

and  let  (3.6)  hold.  Then 


(4.2) 


S„<x»t)  =  l 


( Y+  2wik)x 


"  (Y  +  2itik)n+1 


/  L 


(Y+  2*ik) 


n+1 


An  alternative  derivation  of  (4.2)  uses  the  recursive  relation  (3.5)  and  proceeds  as 
follows.  To  simplify  notations  we  write 
(4.3)  Sft(x)  -  Sn(xjt)  . 

)c 

Observe  that  Sj(x)  is  the  linear  spline  that  interpolates  the  sequence  (t  ),  and  so 

S^(x)  •  1  ♦  (t-1)x  if  0  <  x  <  1 


However  Sj(x)t~x  is  periodic  with  period  1,  because  Sj (x+1 ) t~x~*  =  tSj(x)  : 

S^xjt"*  »  S^xje  Yx.  Let  its  Fourier  series  be  S^xje  Yx  -  J  a^  e2lti,CX.  For  its 

k 

coefficients  we  find  by  an  integration  by  parts  (see  [10,  §5))  that 


% 


(t-ir 

2*t 


1 

(Y+  2wik)2 


and  therefore 
(4.4) 


v*>  •1SS-L 


2  *  (Y+  2*ik)x 

e _  _ 


(Y+  2*ik) 

Since  S ^ ( x)  is  up  to  a  non- vanishing  factor  identical  with  9^(x;t), 
(3.5),  to  perform  the  operation  /^(Mdu  on  S^(x)  a  total  of  n-1 
divide  the  final  result  by  its  value  at  x*0.  However 


it  suffices,  by 
times,  and  to 


(4.5) 


rx+1  (Y+  2wik)u,  _  eY-1  (Y+2wik)x 
Jx  e  du“mKe 


-7- 


Performing  the  operations  as  described  on  (4*4),  we  obtain  the  fraction  (4.2)* 

5.  S  ( x) t )  +  tX  as  n  ♦  •  fc*  Let  us  assume  (3*6)'  so  that  t  is 

non-negative,  and  write  (4.2)  as 


(5.1) 


S  (xjt)t  X  *  l 
n 


1  2»ikx 

e 


.  n+1 


*/L 


1 


-•  ( Y+2*ik)*‘  '  '  -•  (Y+2*ik) 

n+1 


n+1 


We  multiply  each  of  the  two  series  of  this  fraction  by  Y  ;  except  for  the  two  terms 
for  k  -  0,  which  are  *  1  in  both  series,  the  k-th  terms  in  both  series  are  in  absolute 
value  -  |Y/(Y+2wik) |n+1.  Writing 

(5.2)  P  *  log|t j,  hence  Y  *  log  t  -  P  +  ia  , 

we  find  that 

2  2 

|Y/nr+2irik)|2  -  |(p  +  la)|2/|p  +  i(a+2wik)  | 2 - - 


p2+(a+2nk)2 


From  -w  <  a  <  *  we  find  that 

(5.3)  |a  +  2*k |  >  2ir  -  |a|  >  |a| 
and  therefore 

(5.4)  max  lr/(r+2*lk)l  -  (  P - - rj^2  -  «  <  1 

k  pz+(2ir-|a|r 

by  (5.3).  Moreover  6^  >  0,  because  p  and  a  can  not  both  vanish,  as  we  assume  that 
t  *  1. 

Now  it  should  be  clear  that  the  right  side  of  (5.1)  is  *  1  +  and  that  we  may 

write 

S  (x?t) 

(5.5)  - 1  ■  0(6”  )  uniformly  for  x  6  R  . 

tx  t 

Notice  that  the  approximation  (5.5)  deteriorates  as  a  approaches  ±  *  because  6 
approaches  1.  As  the  constant  in  front  of  the  0-term  of  (5.5)  depends  on  t,  but  not 
on  x  and  not  on  n,  we  have  established  the  proposition 
IV.  I£  t  is  not  negative,  then 

(5.6)  lim  S  (xjt)  -  t*  -  1 1 1 X  ei<XX  for  all  real  x  . 

n  1 


-8- 


Can  the  approximation 


6*  Thecorngutation  of ^the  exponential^  unction  f(x)  *  2X. 

(2.13)#  for  t  «  2,  be  used  to  compute  2X  in  view  of  the  convergence  theorem  (5.6)? 
[7,  §11] ,  and  again  in  [10,  §5)  I  stated  that  this  seems  practicable.  Actually  we  find 
that  this  method  does  not  compete  in  accuracy  with  the  modern  approximations  of  2X  in 
[0,1]  by  appropriate  rational  functions  (see  [5]).  However,  we  will  show  that  by 
appropriate  binary  subdivisions  of  [0,1]  our  approach  becomes  competitive. 

We  introduce  the  natural  number  r  and  change  variables  by 


(6.1) 

defining  F(z)  by 

(6.2) 

For  the  base 
(6.3) 

we  have  F(z)  =  tz. 


1 


F(z)  =  2 


z/2 


f  ( x ) 


t  =  21/2 


and  this  we  can  approximate,  in  view  of  (5.6),  by 


(6.4) 
Setting 

(6.5) 


1  /2 

F(z)  -S  (2;  2  /  ) 

n 


z  =  2r  x  =  [2rx]  +  0  »  v  +  0,  (0  <  0  <  1)  , 


where  [•]  has  its  usual  meaning,  and 

(6.6)  V  -  [2rx] 

However 

(6.7) 

and  by  (6.2)  and  (6.4)  we  have 

(6.8) 

We  recall  that  by  (2.13)  we  have 

(6*9)  S  (0?t)  =  P  (0),  (0  <  0  <  1)  , 

n  n,  r  = 

where  Pn^r(u)  is  the  polynomial 

rn,r,u'  *  lu  " 

(6.10) 


S  (z;21/2  )  =  S  (v+0,t)  =  tVS  ( 0;t )  »  2V/2  S  ( 6 ; 2 1 ) 
n  n  n  n 


2*  -  2V/2rs  (0,21/2r) 

n 


P  (u)  »  (un  +  (?)a.(t)un"1  +  •••  *  a  ( t ) }/a  (t) 
n, *  it  n  n 


,  .  (n)  A  .  (n)  n 

1  +  c,  u  +  •  •  •  +  c  u 
1,r  n, r 


In 
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whose  coefficients  are  compute  by  Euler's  algorithm  (2.7). 

*  8 

How  close  does  Sft(8|t)  approximate  t  ?  In  (5.4)  we  have,  by  (6.3)  and  (5.2),  that 
o*0  and  p  *  log | t |  *  (log  2)/2r,  and  so,  dropping  the  term  *  1  in  the  denominator. 


we  have 


and  (5.5)  shows  that 


1  +  (fift)* 


S  (8,*:)  n+1 

-2-§ -  -  1  -  0(±2^i)  . 

t  2»2 


We  conclude  that  the  approximation  (6.8)  will  be  close,  provided  that  either  n  or  r,  or 
both,  are  of  some  size. 

We  need  the  numerical  values  of  the  coefficients  of  the  polynomial  (6.10).  I  am 
indebted  to  Fred  Sauer,  of  the  MRC  Computing  Staff,  for  a  30  place  table  of  the 


coefficients 


cjn*  for  1-1, 


,n i  n  *  3,.. .,20)  r  -  0,1, 


Computing  with  a  hand-held  calculator,  we  record  here  only  their  values  for  r  -  4  and  n 
-  3  to  12  decimal  places: 

(6.14)  c!3*  -  .04332  16979  37,  c*3!  -  .00093  82380  41,  cl3!  -  .00001  38464  49  . 

1,4  2,4  3,4 

These  should  give  in  (6.8)  about  eight  decimal  places  of  accuracy. 

As  an  example  we  compute  f(w)  *  2*.  Here  x  -  z  -  2rx  -  16t  -  50.26548  246,  and 
so  v  -  50,  8  -  .26548  24574.  Since  50/16  -  (25  +  24  +  2)/24  -  2  ♦  1  ♦  2~3  -  3  ♦  1/8, 
(6.8)  becomes 

2*  -  j 50/16  f(6)  „  8.21/8  p  (6)  ' 

3,4 

We  find  from  (6.14)  that  P,  .(8)  -  1.01156  7538,  and  finally 

2*  »  8.82497  7778 

which  is  accurate  to  seven  decimals.  In  this  computation  we  have  approximated  2*  by  the 


interpolating  cubic  spline  Sn(16x;2^^®>  having  24  -  16 


onente  in  the  interval 


[0,1].  In  implementing  (6.8)  it  is  important  to  represent  the  integer  v  -  [2rx]  in 

v/2r 

binary  notation,  as  we  have  done  above,  in  this  way  2  appears  as  a  polynomial  in 


i 


1 


21//2  having  only  coefficient*  0  or  1  •  For  r  ■  4  and  n  ■  9  Sauer '0  table  of  (6.13) 
gives  a  result  accurate  to  21  decimals. 

7.  Cardinal  spline  interpolation.  In  this  last  section  we  sketch  a  solution  of  the 
problem  of  cardinal  spline  Interpolation  [8,  Lecture  4J : 

The  (y^)  being  prescribed,  we  wish  to  find 
(7.1)  S(x)  e  n 


such  that 

(7.2)  S(j)  -  y^  for  all  integers  j 

This  problem  is  trivial  if  n  *  1  and  we  may  assume  n  >  2.  We  further  restrict  the 
discussion  to  the  case  when 

(7*3)  the  data  (y^)  and  the  solution  S(x)  are  of  power  growth  , 

meaning  that  |y^|  *  0|j|Y  as  j  ♦  for  some  Y  >  0,  and  that  |S(x)l  *  0|x|^  as 
x  *  ±«*,  for  some  6  >  0.  We  shall  deal  only  with  the  question  of  uniqueness  of  the 
solution*  because  the  exponential  splines  *^(x;t)  play  a  decisive  role. 

We  need  the  null-space 

(7.4)  S°  -  (s(x)i  S(x)  e  S  ,  S(j)  •  0  for  all  integers  j} 

n  n 


and  state 


Lemma  1 . 
by  the  n  -  1 
(7*5) 
where  the 


_  0  o 

The  linear  subspace  of  5^,  has  the  dimension  n 

exponential  splines 

sv(x)  -  ♦n(x;Ay)  (v  -  , 

are  the  zeros  (1.15)  of  the  polynomial 


1 


and  is  spanned 


11- 1  , 

(7.6)  n  (t)  -  nl  l  Q  ( j+1 )t J  . 

n  0 

That  the  Sy(x),  called  the  eigensplines  of  Sn*  are  elements  of  follows  from 

(1.13)  and  the  functional  equations 

(7.7)  Sy(x+1)  *  *ySy(x),  (V  *  1*...*n-1) 

Indeed*  S^(0)  *  0  implies  that  Sy(j)  *  0  for  all  j.  For  a  proof  of  Lemma  1  see  [8, 
Lecture  4,  § 3 ] • 
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(7.8) 


1 


The  polynomial  II  (t)  being  reciprocal,  we  have 
n 

X  X  -  A  A  -  •  ••  « 
n-1  1  n-2  2 

and  there  is  an  important  distinction  depending  on  the  parity  of  n. 

1.  n  *  2m- 1  is  odd*  From  the  simplicity  of  the  A  it  follows  that 


(7.9) 


A^  <  •••  <  A  <  -1  <  A  <  •••  <  A  <  0  . 

2m- 2  m  m-1  1 


By  Lemma  1  S(x)  e  implies  that 


(7.10) 


n-1 

S(x)  -  l  cy  Sv(x) 
1 


The  inequalities  (7.9)  and  the  behavior  of  the  splines  (7.5)  at  and  -•  implies  that 

(7.10)  is  of  power  growth  if  and  only  if  all  ■  0.  This  is  the  basis  of  our  first 
result 

V.  n  -  2m- 1  is  odd,  and  the  (y ^ )  are  of  power  growth,  then  there  exists  a 

unique  S(x)  e  Sn  of  power  growth  satisfying  (7.2). 

2.  n  -  2m  is^even .  Now  we  find  by  (7.8)  that  the  A^  satisfy 

(7.11) 
so  that 


A„  „  <  •••  <  A  <  A 
2m- 1  m+1  m 


-1  <  A  <  • ••  <  A  <  0 
m-1  1 


(7.12)  S  (x)  -  ♦  (x>-1) 

m  n 

is  one  of  the  eigensplines.  It  satisfies  Sm(x+1)  *  hence  has  period  -  2; 

NOrlund  [6,  Chap.  2,  §16)  attributes  it  to  Hermite  and  Sonine.  The  bounded  Sm(x)  is  a 

counter-example  to  the  proposition  V  for  even  n. 

We  now  abandon  the  class  5  and  consider  the  new  class 

n 

(7.13)  S*  -  (S(x),  Six  ♦  ;)  e  S  )  . 

n  i  n 

This  is  the  class  of  midpoint  cardinal  splines  having  their  knots  at  x  -  j+  — .  Within 

* 

5^  we  have  the  midpoint  exponential  splines  defined  by 

(7.14)  ♦  (xjt)  •  ♦  (x  ♦  if  t) 

n  n  2 

When  are  these  element  of  the  null-space 

(7.15)  Sn°  -  (S(x)i  S(x)  e  Sn»  S ( j )  -  0  for  all  j)  7 
This  depends  on  the  vanishing  of 


*n(°,t)  -  *n<2't)  -  E  t-0n*t<5  -j>  -  2  t  Qn+1<n+1~  j  +i> 


1  tl_V,(1+  I>  ■  t_n  ^  Vi(j  +I)tj 


which  follows  from  (1.3).  This  shows  that  the  eigensplines  of  S  depend  on  the  zeros  of 


the  polynomial 


yt,  -  2ft  n.  I  Qn+1(J  ♦  . 


which  is  called  the  midpoint  Euler-Frobenius  polynomial.  Again  it  is  monic  and  reciprocal 
having  n  simple  and  negative  zeros 


(7.17) 


The  analogue  of  Lemma  1  is 


Mn  <  Vi  <  —  <  **2  <  U,  <  0  . 


Lemma  2.  The  linear  subspace  S  ,  of  S  ,  is  of  dimension  n  and  it  is  spanned  b 


the  n  midpoint  exponential  splines 

(7.18)  sv<*)  *  •n(x,yv)  *  *n<x  +  (V  *  1,...,n)  . 

*0 

It  follows  that  the  general  element  of  S  is 

n 


5 

S(x)  *  l  c  S  (x)  . 


The  role  of  the  parity  of  n  is  now  reversed,  because  again  we  have 

(7.20)  u  y.  *  w  =  •••  »  1 

n  i  n-1  2 

1.  n  *  2ra-1  is^odd.  Now  (7.17)  and  (7.20)  show  that 

(?.21)  “2n-1  <  .«  <  Wm+1  <  »„  *  -1  <  Vi  <  *”  <  **1  <  0  ’ 

*0 

so  that  S  contains  the  periodic,  hence  bounded,  element 
n  — — — — 

(7.22)  S  (x)  *  ♦  (x+  -1) 

m  n  2 

and  there  can  be  no  uniclty  for  the  problem  (7.2)  within  the  subclass  of  S  *  of  power 


2*  n  -  2m  is^even.  Now  (7.17)  becomes 

3)  U2m  <  •••  <  V,  <  -1  <  Wm  <  “•  <  »,  <  0 


and  these  inequalities  allow  us  to  show  that  if  <7.19)  is  of  power  growth,  then  all 


cy  ■  0.  These  results  lead  to  the  analogue  of  proposition  V: 

vi.  n  ■  2a  is  even,  and  <y^)  are  of  power  growth,  then  there  exists  a  unique 

* 

S(x)  e  Sn  of  power  growth  satisfying  (7.2)* 

Concluding^renarjcs .  1.  The  special  case  when  the  data  (y^)  are  bounded,  and  the 

solution  S(x)  is  to  be  bounded,  the  propositions  V  and  VI  were  first  established  by 
Subbotin  [11],  2.  If  in  (7,2)  we  have  y^  ■  (-1)^  for  all  j,  then  the  solution  of  (7.2) 
is  the  function 

(  ♦  (xf-1)/#  (0,-1)  if  n  is  odd  f 
in  n 

(7.24)  E  (x)  -  <  .  . 

I  ♦  <x+  -i-1)/9  (-1-1)  if  n  is  even  • 

>■  n  2  n  2 

It  is  called  the  Buler  spline  of  degree  n,  and  it  is  the  solution  of  the  famous  Landau- 
Kolmogorov  extremum  problem  (see  [9]  for  references). 
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